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We present the results of an investigation of a mixed action approach of overlap valence and 
maximally twisted mass sea quarks. Employing a particular matching condition on the pion mass, 
we analyze the continuum limit scaling of the pion decay constant and the role of chiral zero 
modes of the overlap operator in this process. We employ gauge field configurations generated 
by the European Twisted Mass Collaboration with linear lattice size L ranging from 1.3 to 1.9 fm. 
The continuum limit is taken at a fixed value of L = 1.3 fm, employing three values of the lattice 
spacing and two values of the pion mass constructed from sea quarks only. 
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Lattice size 


L [fm] 




a [fm] 


aji 


m n [MeV] 


#conf 


16 3 x 32 


1.3 


3.9 


0.079 


0.004 


300 


544 


20 3 x 40 


1.3 


4.05 


0.063 


0.003 


300 


300 


24 3 x 48 


1.3 


4.2 


0.051 


0.002 


300 


401 


16 3 x 32 


1.3 


3.9 


0.079 


0.0074 


450 


260 


20 3 x 40 


1.3 


4.05 


0.063 


0.006 


450 


299 


24 3 x 48 


1.3 


4.2 


0.051 


0.005 


450 


137 


20 3 x 40 


1.6 


3.9 


0.079 


0.004 


300 


239 


24 3 x 48 


1.9 


3.9 


0.079 


0.004 


300 


435 



Table 1: We give the simulation parameters and indicate the values of L, a and m n to provide estimates for 
the physical situation. 



1. Introduction 

Mixed action simulations with overlap fermions are considered to be a cost-effective alterna- 
tive to fully dynamical overlap simulations. They allow to profit from the good chiral properties of 
overlap fermions in the valence sector, while in the sea sector a cheaper fermion discretization is 
used to keep the simulation time to a tolerable level. 

We have used several ensembles of gauge field configurations generated by the European 
Twisted Mass Collaboration (ETMC) to investigate the particular setup of overlap valence quarks 
on a maximally twisted mass sea. Our aim was to explore the potential of this approach and study 
the continuum limit behaviour of this mixed action setup. For earlier results regarding this setup, 
we refer to % |, |]. 

2. Simulation setup 

We use ETMC Nf = 2 dynamical maximally twisted mass (MTM) gauge field configurations 
and we refer to refs. [Q, ||] for the details concerning their generation. The parameter values are 
gathered in Tab. |l[ We focus mainly on small- volume ensembles with L pa 1 .3 fm, with pion masses 
corresponding to about 300 and 450 MeV, but we also consider larger volumes to investigate the 
size of finite volume effects in the case of the coarsest lattice spacing. 

In the valence sector, we use overlap fermions defined by: 

Awerlap(O) = \ (l -A^U)" 1 / 2 ) . (2.1) 

with the kernel operator A = 1 + s — aDwiison(O), where s is a parameter which satisfies \s\ < 1 and 
can be used to optimize locality properties and the Wilson-Dirac operator is defined by: 

OwikonM = 2 iMK + V m) -^V* V^) +m, (2.2) 

where m is the bare Wilson quark mass and , V* are the forward and backward covariant deriva- 
tives, respectively. The massive overlap Dirac operator is given by: 

( am q \ - 

^overlap (m q ) = [l + s — J D overlap (0) + m q , (2.3) 
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Figure 1: (left) Decay rate of the overlap operator, (right) Continuum limit scaling of the ratio m K /p. 



where m q is the bare overlap quark mass. 

For each ensemble, following ref. |0] we use the SUMR solver with adaptive precision and 
multi mass capability to compute the propagators for a wide range of valence quark masses m q , 
from quark mass matched to the unitary one to around the mass of the physical strange quark. 
Before applying the overlap operator, we perform one iteration of HYP smearing [Jsj] on the gauge 
field configurations. 

3. Locality of the overlap operator 

In order to achieve the best locality properties of the overlap operator, the parameter s can 
be tuned. It has been shown [Q] that the overlap operator is local for a wide range of simulation 
parameter values, i.e. the norm of the overlap operator falls off exponentially \\D over i ap \\ max {d) 
e~ pd , where d is the taxi-driver distance between lattice points and || . || max is the same norm as used 
in ref. [||]. We have computed the decay rate p for several values of the parameter s, both for the 
overlap operator constructed on HYP smeared and on original, unsmeared configurations. 

We find (the left panel of Fig. [I]) that in the HYP-smeared case the best locality is observed in 
the neighbourhood of the free-field optimal value s = 0. Following this observation, we set s = 
in our computations. We also investigate the continuum limit of the ratio m n j p. At finite lattice 
spacing the condition m % < p must hold [ |l0| ] in order that the interaction can be considered local 
from the point of view of the considered particle. For the pion, the ratio m K j p (right panel of Fig. 
[I]) is well below 1 and thus locality is guaranteed. Moreover, the continuum limit value of m n j p is 
consistent with zero, which is related to the fact that m % had been fixed to m n « 300MeV and has 
a non-zero value in the continuum, while 1 jp vanishes for a = 0. 

4. Continuum limit scaling test of the pion decay constant 

To minimize the unitarity violations present in the mixed action setup, a procedure of matching 
of the quark masses has to be performed. In general, there is a large freedom to choose a matching 
condition. In this proceeding, we will discuss a particular matching condition which will illusttate 
the role of chiral zero modes of the overlap operator especially. Our condition consists of matching 
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Figure 2: The quark mass dependence of the pion mass (left) and the pion decay constant (right). Ensemble: 
j8 = 3.9, L/a = 16, m n « 300 MeV. 
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Figure 3: The continuum limit scaling of the pion decay constant (extracted from the PP correlator) - 
overlap pion at the matching mass and two other reference values of rom^ (left) and MTM pion (right). 



the value of the charged pion mass , i.e. to enforce that the pion built from two valence quarks has 
a mass closest to the one built from two sea quarks. The procedure of matching is illustrated in the 
left panel of Fig. ^, which depicts the quark mass dependence of the pion mass. We extract the 
pion mass from two correlators - the pseudoscalar one (PP) and the difference of the pseudoscalar 
and the scalar ones (PP-SS). Taking the latter, one can avoid any effects from the chiral zero modes 
of the overlap operator, since the zero modes couple equally to the PP and SS correlators. Thus, 
we calculate two matching quark masses for each ensemble. 

It is clear that the effects of zero modes are non-negligible for this ensemble and the two 
definitions of the matching mass lead to different values of the matching quark masses. Hence, we 
will use both definitions and analyze the effects at both matching quark masses. 

The right panel of Fig. || shows the quark mass dependence of the pion decay constant f n , 
where this observable has been extracted again from both the PP and the PP-SS correlator. The two 
curves tend to each other for large values of the quark mass, where the effects of the zero modes 
are small, but for small quark masses (including the matching mass) these effects are rather large. 

In the following, we concentrate on the continuum limit scaling test of the pion decay constant. 

'in the twisted mass formalism, the neutral and charged pions have unequal masses, which is due to 0(a 2 ) isospin 
symmetry breaking. We choose to match the charged pion mass, because of its smaller statistical uncertainty and smaller 
discretization errors. 
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Figure 5: (left) The continuum limit scaling of the difference ro(f% ver ap — fff™) in the case of a heavier 
sea quark mass (m % ss 450 MeV). (right) The safe, hazardous and non-safe regions of parameter space with 
respect to the zero modes effects. 



We will perform it at three values of the pion mass, corresponding to r^m^ ss 0.85, 1.0, 1.5, where 
the smallest value corresponds to the matching mass. Fig. |3| shows the scaling of f n (extracted 
from the PP correlator) for overlap valence quarks (left) and MTM valence quarks (unitary setup; 
right). Clearly, both discretizations show &{a 2 ) leading discretization effects. 

However, when employing the matching condition discussed above, the continuum limit of f % 
in the mixed action setup is different from the unitary one, if f K is extracted from the PP correlator 
(left panel of Fig. [|) at the values of the lattice spacing used here. This can be attributed to the 
effects of the zero modes, since the continuum limit value of f n extracted from the PP-SS correlator 
agrees nicely with the unitary value (Fig. ||, right). We note that the seemingly inconsistent con- 
tinuum limits originate from the particular matching condition used here. Employing alternative 



matching conditions can substantially alter this effect [ 1 1 



The hypothesis about the role of the zero modes can be confirmed by performing an analogous 
scaling test for a heavier value of the sea quark mass (corresponding to a pion mass of about 450 
MeV). In this case, the zero modes effects (which enter the PP correlator as l/ntg and l/m q effects) 
should be reduced. Indeed, the left panel of Fig. || shows that the mixed action and unitary values 
of f n agree in the continuum limit, both in the PP and the PP-SS case. The effects of the zero 
modes are still visible, since f K extracted from PP and from PP-SS slightly differ, but these effects 
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Figure 6: (left) The quark mass dependence of the nucleon and delta mass. Ensemble: j3 = 3.9, L/a= 16, 
a^ = 0.004. (right) The unitary (MTM) and mixed action (overlap) light baryon masses vs. a 2 . 



are strongly suppressed with respect to the m n 300 MeV case. Hence, when using a mixed action 
setup with overlap valence quarks, special attention has to be paid to the effects of the zero modes. 

The investigation of finite volume effects at the coarsest lattice spacing, together with the anal- 
ysis of the sea quark mass dependence at small volume, allowed us to find the range of parameter 
values that ensure that one is safe against the effects of zero modes. This is summarized in the right 
panel of Fig. |5| In terms of m n L the safe region corresponds to m n L > 4 and the non-safe one to 



m K L < 3. For the details of this analysis, we refer to an upcoming publication [11]. 



5. Light baryon masses 

In order to investigate the effects of the mixed action setup in other observables, we have 
computed the light baryon masses in the mixed action and the unitary setup. We have used smeared- 



smeared correlators with the same setup as in Q13 ] . The quark mass dependence of the nucleon and 



delta masses for one ensemble are depicted in Fig. |6| (left). At the matching mass (indicated by 
the vertical bar), the masses in both setups are compatible with each other, within statistical error. 
The situation is similar at a smaller lattice spacing, which is shown in the right panel of Fig. ||. 
This indicates that the zero mode effects are much smaller in the baryon sector, at least in the case 
of light baryon masses. For more details about this analysis, we refer to an upcoming publication 



[12]. 



6. Unitarity violations in the scalar correlator 

The final aspect that we shortly discuss are the unitarity violations in the scalar correlator, 
present in the mixed action setup. As we have already stated, they can be minimized by matching 
the quark masses, but even at the matching mass they can not be entirely eliminated. One effect of 
this kind regards the scalar correlation function, which can obtain an unphysical negative contri- 
bution from one kind of diagrams. This effect in the mixed action setup has been analyzed within 
the framework of Partially Quenched Chiral Perturbation Theory in ref. [|14|]. At large time t, the 
dominant contribution to the scalar correlation function is: 



R 2 „-2M vv t 
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Figure 7: (left) The SS subtr. correlation function at the matching mass . (right) Continuum limit scaling of 
the fitting parameter y. 



If we define y = 2 ^^ 3 a 2 and take periodic boundary conditions in time, we obtain: Css(t) ' =°° 
— y(te~ 2Mvvt + (T — t)e^ 2Mvv ^ T ^y In order to analyze this formula, we use the SS correlator 
with explicitly subtracted zero modes - a field theoretically not clean procedure which we use here 
only to have some estimate of the unitarity violation effect. We then extract 7 for each light-quark, 
small- volume ensemble, by fitting formula ( |6?T| ) to the averaged SS correlator. An example of 
such fit is depicted in Fig. [7] (left). The right panel of this figure shows the extrapolation of 7 to 
the continuum limit. The linear dependence in a 2 and the value consistent with zero for a = is 
compatible with the hypothesis that the unitarity violations in the mixed setup are iff (a 2 ) effects 
and can be described by formula (|Q|). For more details about this test, we again refer to @. 
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